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ABSTRACT 

The discovery of the first gravitationally redshifted spectral line from a neutron star (NS) by Cottam, Paerels 
and Mendez has triggered theoretical studies of the physics of atomic line formation in NS atmospheres. Chang, 
Bildsten and Wasserman showed that the hydrogenic Fe Ha line formed above the photosphere of a bursting 
NS is intrinsically broad. We now include rotational broadening within general relativity and compare the 
resulting profile to that observed during Type I bursts from EXO 0748-676. We show that the fine structure 
splitting of the line precludes a meaningful constraint on the radius. Our fitting of the data show that the line 
forming Fe column is log[()(A^Fe,n=2/cm"^) = 17.9;!;q42 gravitational redshift 1 +z = L345;!;QQQg with 95% 
confidence. We calculate the detectability of this spectral feature for a large range of spins and inclinations 
assuming that the emission comes from the entire surface. We find that at 300 (600) Hz only 10-20% (5-10%) 
of NSs would have spectral features as deep as that seen in EXO 0748-676. 

Subject headings: stars: abundances, surface - stars: neutron - X-rays: binaries, bursts - lines: formation 



L INTRODUCTION 

The observation of gravitationally redshifted atomic ab- 
sorption lines from the surface of the neutron star (NS) in 
EXO 0748-676 (Cottam, Paerels & Mendez 2002; hereafter 
CPM) provides a new constraint on the nuclear equation of 
state. By adding spectra of 28 Type I X-ray bursts observed 
with XMM-Newton, CPM detected an absorption line at 13 A 
during the peak of the burst, which moved to 13.75 A as the 
surface temperature declined. They associated these features 
with hydrogen-like and helium-like Fe « = 2 ^ 3 transitions, 
implying a gravitational redshift of 1 +z « 1 .35. Motivated by 
these observations, Bildsten, Chang & Paerels (2003; here- 
after BCP) explained why Fe can be present above the photo- 
sphere. They also discussed the microphysics of the resonant 
Fe Ha line, showing that it is Stark broadened because of the 
large photospheric densities, p w 0. 1 - 1 g cm"-' (Paerels 1997) 
and confirmed London, Taam & Howard's (1986) conclusion 
that the radiation field determines the Fe ionization balance. 

Chang, Bildsten & Wasserman (2005, hereafter CBW) re- 
cently amplified these points and performed a self-consistent 
radiative transfer and statistical equilibrium calculation. They 
assumed that the Fe Ha line forming region is a thin layer 
above the continuum photosphere, as expected from BCP's 
accretion-spallation scenario. The Ha line is a scattering 
dominated resonance line with a fine structure splitting of 20.7 
eV comparable to the observed width and the scale of rota- 
tional broadening for a 10 km NS rotating at 45 Hz of 24 eV. 
The lack of pulsations from this NS during accretion leads us 
to estimate that B<l-2xlO''G needed for Zeeman splitting 
to play a role in the formation of the Fe Ha line (Loeb 2003). 
CBW calculated the equivalent width (EW) as a function of 
Fe column, A^pe, for varying effective temperatures and pres- 
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sures in the thin line-forming region. To match the observed 
EW required, a total iron column of A^pe ~ 1 - 3 x 10^" cm"^ 
was needed, higher than predicted by BCP by a factor of a 
few. 

In this Letter, we rotationally broaden the energy and 
angle dependent line profiles of CBW and compare them 
to the observations. We explain our method of cal- 
culating the rotational broadening in §2 and show that 
the Schwarzschild-i-Doppler (S-i-D) approximation (Miller & 
Lamb 1998; Braje, Romani, & Rauch 2000) is extremely ac- 
curate in reproducing the line profiles from the full general 
relativistic calculation, even up to 300 Hz. We fit the data 
in §3 and show that the intrinsically broad line profile pro- 
hibits any meaningful constraint on the NS radius if the 44.7 
Hz burst oscillation seen by Villarreal & Strohmayer (2004) 
is the spin frequency. At slow spins or at low inclinations, the 
line profile shape is set by both rotational broadening and fine 
structure splitting, making it absolutely necessary to use real- 
istic line profiles when modeling the lines in EXO 0748-676 
or any other favorably oriented (and hence detectable) rapidly 
rotating NS. We confirm the importance of the intrinsic line 
profile in determining NS parameters by constrasting it with a 
simple line profile. We close §3 by accurately determining the 
redshift and the amount of Fe in the upper atmosphere. We 
conclude in §4 by summarizing our findings and discussing 
the likelihood that lines like that seen in EXO 0748-676 will 
be found in the more rapidly rotating NSs. 

2. ROTATIONAL BROADENING 

We compute the impact of rotational broadening on CBW's 
intrinsic line profiles using two independent methods: a fully 
relativistic method that employs ray-tracing in the spacetime 
of a rotating NS and an approximate method which retains the 
most important relativistic effects to lowest order in v/c. 

In the fully relativistic approach, the spacetime metric 
for rotating NSs with various spin frequencies, masses and 
equations of state are computed numerically as described in 
Cadeau, Leahy & Morsink (2005a). Once the (non-spherical) 
surface of the NS has been located, geodesies connecting the 
surface to the observer are computed. We discretize the sur- 
face into small patches and calculate the zenith angle a be- 
tween the local normal and initial photon direction in the lo- 
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cally comoving frame for each patch. We also calculate the 
solid angle dfl subtended by the patch as seen by the observer 
and the redshift z* ■ The redshift, z* includes both the gravita- 
tional redshift and Doppler-like terms resulting from rotation. 
Since the intensity, /, is both a function of energy and cos a, 
the observed flux due to the patch is 
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The observed spectrum is found by summing over all visible 
patches. The technical details of the fully relativistic compu- 
tational method are presented by Cadeau et al. (2005b). 

We also use the simpler S+D approximation (Miller & 
Lamb 1998; Braje, Romani, & Rauch 2000; also see Cadeau 
et al 2005a and Poutanen & Gierlinski 2003) to indepen- 
dently compute the spectrum. In the S+D approximation 
the metric is approximated by the non-rotating Schwarzschild 
metric, and appropriate Doppler factors correct for the ro- 
tation. To lowest order in v/c, the redshift factor reduces 
to (1 + z*) = (l+z)(l-(5) where S = /Ssin^sinasin/, /3 = 
27r7?i/s/(c\/l - 2GM/(Rc^)) is the scale of rotational broad- 
ening, is the observed spin of the NS, and z = (1 - 
2GM /(Rc^))~^l^ — 1. The angle C, is the azimuthal angle of 
the emitting patch about the vector pointing from the center 
of the star towards the observer This angle is zero if the patch 
is in the plane defined by the spin axis and this vector 

A further simplification is the use of the Beloborodov 
(2002) approximation, 1 -cosa = (1 +z)"-(l -cos?/;), where 
ij] is the bending angle. This allows the total observed flux to 
be written as the integral over all possible initial photon zenith 
and azimuthal angles, 

-if^Vf^^f^ ^ I (1^(1 +z.), cos a) 
d) J J ii^^ ' 
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where R is the NS radius and d is the distance. The results for 
emission from the entire NS surface with CB W's intrinsic line 
profile of the fully relativistic calculation and the S+D approx- 
imation are shown in Figure [0 They are in excellent agree- 
ment and the tiny deviations are mainly due to our neglect of 
relativistic aberation. The linear change to the metric in v/c 
due to frame dragging (Hartle 1967, Hartle & Thorne 1968) 
suggests observable effect on the line at spin frequencies near 
300 Hz as /3 ~ 0.1. However, the frame-dragging term enters 
into the line broadening calculation as a second order effect in 
v/c, preserving the validity of the simple S+D approximation. 
Only at very high spin frequencies does frame dragging have 
the potential of producing observable effects (Bhattacharyya 
et al 2004). 

Of greater relevance is the relative size of the fine struc- 
ture splitting and the scale of rotational broadening. At 45 
Hz, the scale of rotational broadening is 24 eV for a 10 km 
NS, which is similar to the fine structure splitting of 20.7 eV 
for the Fe Ha line. At slower and faster rotation rates, fine 
structure splitting dominates rotational broadening and vice 
versa. Hence the line profile can change dramatically near 45 
Hz. We illustrate this effect in Figure^where we plot the line 
profile for z^s sin / = 20 Hz (long-short dashed line), which is 
dominated by fine structure splitting. 

3. FITS TO EXO 0748-676 AND IMPLICATIONS 

The simplicity and accuracy of the S+D approximation al- 
lows for a rapid and detailed exploration of parameter space. 
We assume emission from the entire surface. Emission con- 
centrated near the spin axis would also give narrow line pro- 
files (Bhattacharyya, Miller & Lamb 2004), but we do not 




Fig. 1 . — Line profiles from an exact GR calculation (thick lines) and the 
S+D approximation (thin Hnes). We plot four different models with Us sin i = 
45 Hz (solid lines), j^ssin! = 100 Hz (short dashed lines), j^ssin; = 200 Hz 
(dotted lines), and i/f, = 300 Hz (long dashed lines). The intrinsic line profile 
(dot-dashed line) is the same for these models. We also plot the line profile 
for i/f, sin; = 20 Hz (long-short dashed line) to illustrate the fine profile as it 
changes from being dominated by fine structure splitting to being dominated 
by rotational broadening. 



consider this possibility. We fit our models to the co-added 
spectra and continuum model from CPM. We focus on a nar- 
row range between 12A and 14A and plot these points in Fig- 
ure|2] The 13.25 A feature in the continuum fit is due to a Ne 
IX resonant transition in the circumstellar model. Taking the 
continuum model as the background flux, we overlay some 
models for = 0, 45, 100 and 300 Hz for R=10 km. To set 
the scale of Stark broadening, we take a background density 
of Hp = 10^^ cm"^, which is appropriate for a uniform Fe abun- 
dance above the continuum photosphere (CBW). From these 
models, large spins i.e. i^ssin/ = 300 Hz are ruled out. Sec- 
ondly, the intrinsically broad line profile which is similar to 
the width of the observed line, prohibits a meaningful limit 
on R for = 44.7 Hz. 

We also performed a detailed exploration of parameter 
space. Using the data between 12.6 A and 13.3 A, we overlay 
our line profiles on the continuum model to test the goodness 
of fit. Varying t/, sin /(/?/ 10km), 1h-z, and A^Fe,n=2, the best fit 
gives a = 11.01 for 17 degrees of freedom. By itself, the 
continuum model gives = 24.9. The best fit parameters 
arej/sSin/(7;/10km) = 32±19Hz, l+z= 1.345±0.004, and 
logio(A?Fe,n=2/cm-2)= 17.9 ±0.2. The errors indicated are one 
sigma values derived from the covariance matrix. 

To develop a more detailed understanding of parameter 
space, we marginalize our three parameter joint probability 
distribution function (PDF) over two parameters at a time and 
plot the resulting one dimensional PDF for each parameter in 
Figure|3l The 68%, 95%, and 99% central confidence inter- 
vals for the intrinsic line profiles for sin /(/?/ 10 km) (Rsini 
for j/s = 44.7Hz) are ll-52Hz (0.25-llkm), 4-84Hz 
(0.09- 19km), and 2- 122 Hz (< 27km) respectively. For the 
best measured value of 1 + z = 1.345, a 1.4 Mq NS with R = 
9.2 km falls within these limits. 

The PDF for i^, sin /(7?/ 10km) is asymmetric toward zero 
and hence prohibits setting a meaningful lower limit due to 
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Fig. 2. — Line profiles atR= 10 km for several different models compared 
to the data. We plot best fitting line profiles for = Hz (dotted line), 44.7 
Hz (short-dashed line), 100 Hz (long-dashed line), and 300 Hz (dot-dashed 
line) for if sin; = 10 km. We also show the continuum model of CPM (solid 
line) and the data (one sigma error bars). 



the instrinsic width of the line set by its fine structure. Had 
we ignored the instrinic shape of the line and had adopted a 
black line profile, which is a narrow line with zero flux about 
line center and fixed EW = 0.1 A (dashed Une), we would set 
misleading constraints on the neutron star parameters. Using 
this black line profile, we can constrain ly sin i(R /10km) = 47- 
102 Hz which for a 44.7 Hz rotator gives Rsini = 10.4-23 
km. These misleading constraints highlight the importance of 
using the intrinsic Une profile to fit for NS parameters. The 
limitations of simple line profiles have also been recognized 
previously by Villarreal & Strohmayer (2004). 

The 68%, 95% and 99% central confidence intervals of 
logio(A'Fe,n=2/cm-2) are 17.69- 18.05, 17.48- 18.17, and 
17.34- 18.26 respectively. The corresponding A^pe's are 
1.4-3 X lO^^cm-^ 0.9-4.0 X lO^^cm'^ and 0.8-4.8 x 
lO^^cm"^. At 68% confidence, these values are within a 
factor of 3-6 times the solar metallicity photospheric value. 



A^Fe ~ 5 X 10'^ cm"^ (CBW) and within a factor of 5-10 times 
the accretion-spallation value of A^Fe ~ 3.4 x lO'^cm"^ (BCP). 
At 99% confidence, we may require a colunm a factor of up 
to 10 times larger. Our estimate from this likelihood analy- 
sis agrees with CBW's simple EW estimate. In addition the 
NLTE value of A'fb is likely larger (London, Taam & Howard 
1986) and may demand either supersolar metalhcity or radia- 
tive levitation (CBW). 

The 68%, 95% and 99% central confidence intervals of 
1+z are 1.342-1.348, 1.337-1.350, and 1.334-1.353. We 
have also calculated the corresponding confidence intervals 
for R/M which gives R = 6.592-6.67O(M/M0)km, 6.565 - 
6.722(M/MQ)km, and 6.542 -6.766 (M/M©) km. Since we 
fit a rotationaUy broadened profile to the data, our measure- 
ment of 1+z is immune to beaming effects which Ozel & 
Psaltis (2003) (also see Bhattacharyya, Miller & Lamb 2004) 
point out can skew the measurement of the redshift if the line 
energy is taken from flux minimum (or combinations of min- 
ima). 

4. SUMMARY AND CONCLUSIONS 
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Fig. 3. — Marginalized probability distribution functions for 1+z, 
i^s sin!(R/10km) and ^^-5.11=2 normalized so that their peak is unity. For 
the marginalized PDF of i^s sini(R/10km), we also display on the top axis 
the i?sini for = 44.7 Hz. In addition, we plot the marginalized PDF of 
j^s sin ('(/?/ 10 km) for a black line profile with fixed EW of 0.1 A (dashed 
line). We also show the 95% (dotted line) central confidence intervals for 
each parameter. 



Using the intrinsic Fe Ha line profiles calculated by CBW, 
we calculate rotational broadened line profiles to compare to 
file observed Unes in EXO 0748-676 (CPM). We show fliat the 
Une profiles generated by the simple S+D approximation ac- 
curately reproduces the ones from a fully relativistic ray trac- 
ing method described briefly in § 2. 

Assuming that the burst emission comes from the en- 
tire surface of the NS, we constrain the Fe column in the 
line forming region (log[QA^Fe,n=2/cm"^ = 17.9;!;q42). '^^^^ 
find that the redshift is 1 + z = 1.345!jJ;5][J^. The NS spin is 
i^s sin /(/?/ 10km) = 32!^^ Hz for a fiducial NS radius of 10 
km. These errors denote 95% confidence intervals. Fixing 
Us = 44.7 Hz (ViUarreal & Strohmayer 2004), we find that the 
radius is effectively unconstrained by the observed line due 
to the intrinsic width set by fine structure splitting, an effect 
missed by simple line profiles. 

Other transitions would simultaneously confirm the redshift 
measurement and may provide meaningful constraints on the 
radius. CBW calculated the line profiles for the associate Lya 
transition and Pa transition. The EWs of the Fe Lya and Pa 
lines were « 15-25eV and « 7 - 1 e V, respectively. The Pa 
transition is obscured by interstellar absorption. CBW points 
out that the Lya line profile is dominated by rotational broad- 
ening. Hence a sufficiently deep observation of the Lya tran- 
sition would provide a good constraint on the radius of EXO 
0748-676. The helium-like Fe n = 2 ^ 3 may provide some 
additional constraints on EXO 0748-676, but the intrinsic Une 
profile is not known. This is critical in modeUng this feature, 
but is beyond the scope of this work. Spectral modeling of 
the continuum emission during Type 1 x-ray burst has also 
been utilized to determine NS parameters (Shaposhnikov & 
Titarchuk 2004; Shaposhnikov, Titarchuk, & Haberl 2003 and 
references therein). Redshift measurements and spectral mod- 
eling utilized in combination would set very tight constraints 
on NS parameters. 

Though EXO 0748-676 is rotating slowly among LMXBs, 
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Fig. 4. — Fraction of Sky for which the line minimum is below 10% (thin 
lines) and 30% (thick lines) below the continuum as a function of spin. We 
fix the Fe Ha line to the observed EW from EXO 0748-676 and plot the line 
for NS Rao = 10 (solid lines), 15 (dotted lines) and 20 km (dashed lines). 



more rapid rotation does not necessarily make spectral lines 
undetectable. If a rapidly rotating NS is viewed more face 
on, the spectral lines will be narrower Hence, more rapid 
rotators may just constitute a smaller population of NSs with 
detectable features. We illustrate this in Figure |4] where we 
plot the percentage of sky denoted by H/Att over which the 
line minimum exceeds 10% and 30% below the continuum. 
We assume the same intrinsic line profile as the observed Ha 
Une on EXO 0748-676 for a NS Roo of 10, 15 and 20 km. 
We use Roo instead of R because Roc is a convenient way of 
absorbing the factor of 1 -I- z in the scale of broadening, /3. 
These lines follow a very simple trend, which we now derive. 
Up to a certain rotation rate, these lines are always detectable. 
For sufficiently large rotation rates, the FWHM of the line is 
oc i/ssin/. Since the EW is conserved, the line minimum is oc 
(z^s sin /)"' . Therefore, maximum inclination for line minima 
greater than a certain strength is sin /max oc v'^ which limits 



the fraction of sky over which this line would be detectable. 
Hence, the larger effect of rotational broadening with higher 
inclination broadens the line and weakens line minima unless 
emission is restricted to the rotational pole (Ozel & Psaltis 
2003; Bhattacharyya, Miller, & Lamb 2004). 

For a fiducial value of Roo = 10- 15 km, 10-20% of NS 
with spins of 300 Hz would have spectral features similar to 
EXO 0748-676 where the line minimum is about 30% below 
the continuum. At 600 Hz, the percentage is reduced by a 
factor of two. For these favorably spinning and oriented NSs 
(line minimum « 30%), the line profile is set by both rota- 
tional broadening and fine structure splitting. Therefore, burst 
studies of NSs with spins near 300 Hz may yield at least an- 
other example of atomic spectral lines from NSs. For these 
most favorable systems, the intrinsic line profile is absolutely 
necessary in order to model these lines and perform mean- 
ingful measurements from the observations. Of the presently 
known LMXBs with X-ray bursts, a very promising target is 
4U 1728-34 which spins at 363 Hz (Galloway et al 2003) and 
exhibits bursts at a fairly regular rate. A recent model of the 
NS-accretion disk geometry (Shaposhnikov et al. 2003) esti- 
mates / to be around 50°. At this inclination, the line mini- 
mum is 7% below the continuum. For smaller inclinations of 
i = 30° and 10°, the line minimum would be 10% and 26% 
respectively. 
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